The quantum theory of a quadratic gravity action 
for heterotic strings 



Simon Davis * and Hugh Luckock 



School of Mathematics and Statistics 
University of Sydney 
NSW 2006, Austraha 



Abstract. The wave function for the quadratic gravity theory derived from the heterotic string 

-* 

effective action is deduced to first order in ^-3- by solving a perturbed second-order Wheeler- 
DeWitt equation, assuming that the potential is slowly varying with respect to Predictions for 
inflation based on the solution to the second-order Wheeler-DeWitt equation continue to hold for 
this higher-order theory. It is shown how formal expressions for the average paths in minisuperspace 
{(a(t)), (<I>(t))} for this theory can be used to determine the shifts from the classical solutions aci{t) 
and ^ci{t)i which occur only at third order in the expansion of the functional integrals representing 
the expectation values. 
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Given a dimensionless metric and scalar field, the form of a sigma model coupled to gravity [1], 
renormalizable in the generalized sense, is 



6i(</>)(n0)2 + 62(0)(v^0)(v^</>)n0 + 63(0)[(v^<^)(v^0)]' 

+ 64(.^)(V^0)(V''<^) + hicP) + ci((/.)i?(V..^)(V'^</.) 

+ a2{(p)R^,yR'^'^ + a3{(f>)R'^ + ai{(t))R + sur f ace terms 



(1) 



where [64(0)] = 2, [65(0)] = 4 and [a4((/))] = 2, with aj(^) 7^ for at least one i G {1,2,3} and a 
renormalizable potential term b5{(f)). Setting 



\ b,{4>) = -\v{4^) 
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produces an action describing an exponential coupling of a scalar field to quadratic curvature terms. 
Defining the dilaton field to have dimension 1, $ = ^, and V{^) = -^V{(f)), the quadratic gravity 
theory 



R + ^-{B^f + X (^M-«Ai?^""^ - 4i?^,i?'*'' + R') - 



^2" '2' ^' ' 4^1 



(3) 



is equivalent to the one-loop hcterotic string eff'ective action with the coefficient of the RR term 
set to zero. A non-minimal coupling ^R(^^ term, which has been found useful for generating open 
infiationary universe models [2], can be added to this action by setting the coefficient a4(^) equal 
to ^(1 + ^0'). 

A conformal transformation of the R? term produces a Ricci scalar plus an extra scalar field, 
and any factor multiplying R can be eliminated through another conformal transformation [3- 
5]. There is also a Legendre transformation of the R'^^R^^ term to a Ricci scalar together with 
extra tensor modes [6-8]. The scalar and tensor modes contain an intricate dependence on the 
derivatives of the metric and Ricci tensors, which complicates the definition of the corresponding 
momenta, so that there is no simplification in the canonical quantization procedure by applying a 
conformal or Legendre transformation to the action. The R^^p" R^^pa can be rewritten as a linear 
combination of C^^'^f^Cn^pa and terms containing R^^, and R. The renormalizability of the action 
with C^^^'"^Cfj,vprj has been established [9]. Its efi!'ect on unitarity must be counterbalanced by that 
of the other tensor modes. 



For an action with a polynomial coupling between the scalar field and the quadratic curva- 
ture term ai ((/>), 02 (^) and az{4') will be shifted to di{4>), 0,2(4') 03(0). The coefficients 



in a truncated Taylor series representation of ^-2-, Yln=o 



1 (-1)" 

1,! 



(iVo-hl) 



(-1) 



ATn-l-l e" 



1 



A < 1, would be shifted so that the couplings of the renormalized theory are ^n=Q + 

{No\+iy. ^i^°~^^\^' ^) ^ ^' ^ "^^^ remainder term can be made arbitrarily small and the series 

converges as A'^^ — > oo unless Ci^n has a factorial dependence on n, so that the new action with 
coefficients Oj($), i = 1,2, 3, will be defined and renormalizable in the generalized sense. 

After imposing the restriction to the minisuperspace of Priedmann-Robertson- Walker metrics, 
and adding a boundary term to the action to eliminate terms containing a, the one-dimensional 
action is 



dt 



_<j) 

94 



(4) 



where a{t) is the scale factor of the Priedmann-Robertson- Walker universe, which is open, flat or 
closed if = —1,0 or 1. 

Given the conjugate momenta 



— <3> 

Pa = -12aa + e^^a^ + K) 

94 
94 



(5) 



the Wheeler-De Witt equation to first order in ^p- [10] is 

_<j) 

= (Ho + ^Hi)^ = 

94 

" 24 da a da 2a^ d<l>^ ^ ' 

I fK 1 \ d 1 ^2 I d^ 

Hi = ^ -r - + 
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a4 V 4 576a4 J da 576a^ da^ 1728a6 da^ 

1 f7K 35 \ d 1 d^ 1 d^ 1 d"" 

^ \ ~T ~ 576a4 ) d^ ^ 24a^dad^ ~ MaJ da^d^ ^ SMaPdad^' 

Given that = + -^^i, when |y"^y(<I>)| < 1 and \V{^)\ < 1 in Planck units, so that 
the derivatives of the quantum cosmological wave function ^' with respect to $ are negligible, the 
equation becomes ifo^i ~ —Hi'^q, to first order in ^t-. It may be noted that Hq ^t-^*! would 

.94 V 94 / 

contain a term of the form tt-ti^ I ^-^^i I giving rise to a contribution ^-s-ttt^'i — + 

2a^ ^ip~ "dW' t° differential equation at first order in the expansion parameter. However, if this 

term is included, then the solution to the standard Wheeler-DeWitt equation with the ^2 operator 
and each of the <I> derivative terms in Hi would have to be used. When the $ derivative terms 
are discarded from at the beginning of the computation, it is sufficient to consider the differential 
equation for \E'i without the additional term. 
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The correction to the wave function is then given by ^t-\E'i where 

9i 



Ci*oi + C'2VE'o2 - *02 / -^oi^^ada + 24^'oi / vE'02— da (7) 



^'1 

with the Wronskian defined to be 



d d 

W = *oi3-*02 - *02 3-*01- (8) 
da da 



Consistency with symmetries of the theory depends on the choice of boundary condition, which 
also determines the feasibihlty of obtaining an inflationary cosmology. Upon consideration of the 
N=l supergravity theory restricted to the minisuperspace of Bianchi IX metrics, for example, the 
requirement of homogeneity implies a Lie derivative condition on the spinor fields, which defines 
a no-boundary ground state [11]. While the no boundary wave function is defined to be regular 
in the limit a, 0, this property only holds for the tunneling wave function when the operator 
ordering parameter p is less than one. 

The no-boundary wave function is 

^ONB = 7 2T [y) 

when = — 1 or 1. While the no-boundary wave function is defined by a path integral over 
compact four- manifolds, leading to the conventional choice K = 1, the other values of K are 
possible if the range of coordinates in the fiat or hyperbolic sections is finite. While the probability 
amplitude defined by the no-boundary wave function with a positive coefficient in the exponential 

prefactor does not directly imply the existence of an inflationary universe with the appropriate 
e-folding factor, it may be noted that a negative coefficient can be obtained by the other choice of 
sign of {—zq)^ in the asymptotic expansion of Ai{—ZQ), Zq = z{a = 0) = —K ^ as V{^) — > 0, 

12 ^ 

so that as zo —>■ —00, the normalization factor tends to j^{—zo)~^e^3^~^o)^ if K = 1. Even if 
the coefficient is positive initially, the value of V($) would be driven to zero, so that the change 
in sign can be obtained by analytic continuation in the variable V. The change in sign of the 

exponential prefactor does not affect the regularity of the wave function in the a{t) — >■ limit. 

When K = —1, Ai{—zo) -^{zo)~^sin 
tends to nTT, n integer. 
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, SO that the wave function diverges ^ + f 



When i^^ = 0, it is not necessary to rescale a{t) and V($) to obtain the standard form for the 
second-order Whcclcr-DcWitt equation. Without a term proportional to K, one definition of z 
could be {2V)~^ (a'^V) = a'^V^ , but then z{a = 0) would vanish. If the normalization factor is 
chosen to be Ai(—Zc), where Zc = z{ac), it can be shown that the probability distribution is peaked 
at F = if ttc <C 1. The change in sign of z across the V($) = boundary leads to different 
asymptotics for the wave function. If the positive sign is chosen for the exponent in the prefactor 
eT|(-2:c)2 _ QT^al\v\2 ^ jjjf^g^^jQ^ again can be obtained. 
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Given the two independent solutions of the homogeneous second-order differential equation 
i^o*o = 

^'oi = Ai{-z) ^-02 = Bi{-z) 

the Wronskian [14] is 



V \ 6K 



W = Ai(-z)^Bi(-z) - Bi(-z)^Ai(-z) 
da da 

dz 1 aV^ 2 
da Stt ^ ^ 



dz 
da 



Ai(-z)^Bi(-z) - Bi(-z))^Ai(-z) 
dz dz 



and 



(11) 



*i = CiAi(-z) + C2Bi{-z 
727r 
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^ ^ 36 V 36 / ^ ^ 



KAi'(-z) + ^ Ai"'(-z) 
^ ^ 36 V36y ^ ' 



Ai{-z) 



Bi{-z). 



(12) 



Prom Airy's differential equation, it follows that A"'{—z) can be replaced by zAi' (—z) + Ai{—z) 
in the integral, giving 



*i = CiAi{-z) + C2Bi{-z) 



+ 
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Imposing the Hartle-Hawking boundary condition on the corrected wave function implies that 
it must have the same form as the standard wave function, so that the coefficient of Bi(—z) should 
vanish. 

— / — \ ~^ 

After changing the integration variable, ^ = ^1 + (^) ^ -^j dz, the wave function 

^1 can be obtained by evaluating two integrals of the form J ■jj^^^Ai'{—z)Ai{—z){l + k'z) and 



The integral / dz a{z) S^{(j){z))- S^{i^{z)) has the form [A{z)Si,{(j){z))+BS^+i{(j){z))] S^{'4){z)) 
+[C(z)5^((/)(z)) + D{z)S^,+i{4>{z))\S^+i{i^{z)) when Sf,{z) is a cyUnder function [15] if 



a{z) 






B'{z) + 
C'{z) + 
D'{z) - 



4>{z) 



(^^ + 1) 



(/^+i)§4+(-+i) 



+ Di,\z) - A^'{z) 
C{z) + D<^{z) - Ai}\z) 
"^"'^^^^ Z)(z) - B^\z) - C(j)'{z) 



(14) 



This set of coupled differential equations can be reduced to the 3x3 system 



\C{z) 



where 
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V -'P'iz) 

a{z) 
b{z) 
c{z) 



'A{z) 

+ M I 



cP'iz) 

-(l + 2^)0 + ^</,(.)</>'(.) 
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(15) 
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The solution to this system of differential equations is 

5(z) 1 = expl- M{z')dz' ]■{ r b{z')dz' 



C{z) 



J c{z')dz' 

/ fz \ ( r a{z')dz' 

T-^Texp - / M{z')dz' ■ T-^T /" b{z')dz' 

^ ^ \ f c{z')dz' 



(17) 



If T is the matrix which diagonalizes M{z')dz' ^ then 
Tei-r M{z')dz')r^-i ^ exp(-T f M{z'dz')T-^ 



e-^i 





(18) 



for some set of eigenvalues Ai,A2, A3, and 



lA{z) ' 
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//" a{z')dz' 

-^^ I r biz')dz' 

e-^'' / \ J" c{z')dz' 
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(19) 



Since Ai{z) = (fzi) - Ji 



and Ai'{z) = — |z 



-§ (i 
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dz !i^±KllAi{z)Ai'{z) 



{1-kzy 



-II 



I 



dz 



(1 - kz)' 



:[At{z)]' 



(1 + c'z) , 



dz 



I 1 -Z2 / 2 -Z2 

-3 V 3 / -3 \ 3 



/2 3\ /2 3 

3 I 3 y -3 I 3 



, /2 3\ , /2 3\ , /2 3\ /2 3 

— / 1 -Z2 /2 -Z2 + /l -Z2 /2 -2;2 

"3V3 / 3V3 / 3V3 y^vs 



(1 - /cz)2 



with 0(2;) = |2;2 . 

Substituting tliis function of z into the matrix M[z) and integrating gives 
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The eigenvalues are roots of the cubic equation 
+ 3ln z X'^ + 



-z^ln z - -(3/i^ + i/^ + it + 1/ - l)ln^z 
3 4 



27 



+ -(At + i/ + l)2;^/n 2; 



0. 



(21) 



(22) 



Defining the coefficients a, /3, 7 by using the standard form of the cubic equation X^+aX'^+pX+'j 
0, it follows that 



p = (s/j? + u^ + ij, + u + ^ In^z + ^z^ln z 

9 / 2\ 1 

- G/jiu'^ + 3/x^ - - - At - - Q + o (a* + ^ + 2)2;^Zn2-2;. 



(23) 



When /Lt 



I and 1/ 



2 

'3' 



729 81 216 



1728 



In^z 



(23) 



and Q = when z = 1 01 -^w^ + ^w'^ + — = with w = -j^. Since this cubic equation 
has a single positive root at to = 0.08516, Q > for all z > 1. If to < -2.12470, then Q < 0; if 
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-2.12469 < u; < -0.29559, then Q > 0; and if -0.29558 < u; < 0, then Q < 0. The eigenvalues 
are 



Ai 



+ 



l:z^ln'^z + 
6 



--z^ln^z 
6 



^z^lnh + ^z^ln'^z + ^zHn'z - 
729 81 216 1728 



^z^ln^z + -zVz + ^zVz - 

729 81 216 1728 



1 1 
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iTfz 



--z^ln^z + 
6 



^z'^lrv'z + ^zHn^z + ^zHn^'z 

729 81 216 1728 
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— z^ln^z + —z^ln^z + — 2;^Zn^2; —In^ 
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(25) 
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729 81 216 1728 
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—z^ln^z + —z^ln^z + —z^ln^z - — — 
729 81 216 1728 



Similarly, the following expressions are obtained for p, q and Q for the other values of /Lt and 
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-^z^ln^z - -^z^ln^z 
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-z^ln z — — In^z 
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——z In z — ——z In z 
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(26) 
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The transformation matrix T is equal to 



/ ^"11 


V2I 


V3I 




V22 


V32 




V23 


V33 


/ Vil 
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Vi2 




Mi 


\Vi3 





/ 2 ,4 



9 jj. — f ' 



\ (27) 
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^z■' + [\i + ^{l+2^l)ln z-l-^z'^\[\i-l{n+u)ln z] 



where Hi is a normahzation factor, so that 

A{z) = {e-''-vl^v^^ + e-^-vl^v^2 + e-^'vl^v^3) J a{z')dz' 
+ (e"^'vn?^2i + e-^^'y^2^22 + e-^^vl^V23) j b{z')dz' 
+ {e-^'vl^V3i + e-^'vl2V32 + e-^'vl^vsi) J c{z')dz' 

B{z) = (e-^^w^iVii +e-^^v^2'yi2 + e"^^V23«i3) J a{z')dz' 
+ {e-^'v;^V2i + e-^^t;^2^^22 + e-^^v;r^V23) J b{z')dz' 
+ {e-^'v*2^V3i + e-^^z;;2^^32 + e-^^v^^vsi) J c{z')dz' 

C{z) = ie-^'v*s^vn+e-^^vl2vi2 + e-^'vlsvi3) J a{z')dz' 
+ {e-^'v;^V2i + e-^^vl^V22 + e~^'vl^V23) j h{z')dz' 
+ {e-^^vl^V3i + e-^^vl^V32 + e-^^vl^V3i) j c{z')dz'. 
For the first integral in equation (20) 



(28) 



J a{z')dz' 



, 3-2kz 
kz— ; smh 



1-kz 

5 - 4:kz 



kz 



1-kz 



+ X[Vkz ^ , 
ki\ i-kz 



sink 



kz 



1-kz 



2 , / (l-kz)2\ 2 k' 3 
+ —In 7=^ + -j^z^ 



k' ((1-kz)^ 
A; 2 \ 1 + ykz 



c{z')dz' = 



9 
z^ 



1 9 

— z^. 



81 /i- 1/ 

(29) 

The equality of the indices n, v in the second integral gives rise to a singularity in the matrix 
M{z')dz', which can avoided by using the identity Ii,-i{z) — 1^+1(2;) = ^/^(z), to obtain an 



integral containing the products of modified Bessel functions different indices. The correction to 
the standard wave function can be obtained through the substitution z ^ —z in the integrals (20). 



The contribution of the additional integral in equation (12) to the wave function includes Ai{—z) 
so that the functional dependence of the normalization factor Ai{—zo) + Ci ^^Ai{—zo) + ^^/(^o) 

is divisible by Ai{—zo). Given that the normalization factor contains Ai (^K (^) , the conclusions 
concerning the feasibility of predicting an inflationary universe through the quantum cosmological 
wave function are unaltered by the addition of the higher-order curvature terms. 



In the Planck era, the higher-order curvature terms in the perturbative expansion of the string 

effective action have approximately the same magnitude as the Ricci scalar, and similarly it is 

— *i> 

inappropriate to use a truncated form of a series expansion in of the Wheeler-De Witt equation. 
Instead, a closed-form sixth order differential equation, which can be obtained by including the 
conjugate momentum to derivative of the scale factor Pa, and then using the Ostrogradski method 
to define the Hamiltonian for the higher-derivative theory, can be used to define the quantum 
cosmological wave function in the initial era. Given the one-dimensional action 

and the conjugate momenta 

P. 



—<s 



a{a^ + K) 



(30) 



Pn = 



P^ 



dL 

da 
dL 
da 
dL 



dL 

da 



9l 



6 + 
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{a' + K) 



(31) 



a^^ 



the Hamiltonian is 
H = 



Pah + Pa'd + - L 

— $ 1 

-Qa{o? + K) + %^^{a'^ + K)a + -a^$ 
94 2 



3^2 



(32) 



1 



giP^'e^a'^Pa + ^ 



Pi + 



KP' 



and the pseudo-diflferential equation 
equation * 

, d^^ 2 ^"'^^ A 



can be transformed into the partial diflferential 



Ke 



-2* 



\da?d^ ~ dc^d^J 



4 4 
a 54 



3^* 4 



+ ^94e 



da 



da'' 



I 2 -$ 
+ 0^46 



dad^ 



d^ 



+ 



d 



+ 



4a6 



d^ 



93^ 



(33) 



* This equation differs from the equation (11) in reference [14] by a derivative with respect to 
Imposing an additional constraint on the wave function, the derivative term can be eliminated 
through the addition of an extra term in the Lagrangian. 
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By including a potential term in the Lagrangian and discarding terms containing derivatives of 
^ and V{^) with respect to the following equation is obtained 



?0 + "^^^^ + '^''"^S + [^^'^^ + - ^^^'"'^^^^^^ - ^aWvim = (34) 

The corrected wave function in the inflationary epoch can be matched with the solution to 
a third-order partial differential equation in the initial era along a boundary which also must be 
determined by setting the derivatives up to second order in a to be equal. 

While the wave function *(a, has been used to establish whether the curvature-dependence 
defined by the potential favours inflation, to determine the most probable path in minisuperspace 
{a{t), it is preferable to consider the partition function 

Z = J e-^['^W'*Wl d[a{t)]d[^it)] (35) 



which is extremized at the classical solutions 

6a 



ad. 



6Z_ 



0. (36) 



The expectation values {a{t)) and ($(i)) based on the action / are 



(a(t)) 



J n{i)e. -'["(''■*('':r/[«(/)]r/[<I)(0] 



Sa{t)- 



■(5a(t))= 



J ^(t) e-^[«W'^Wld[a(t)]d[$(0] 
/ e-f[«W'*W] d[a(t)]c/[#(t)] 



d[a{t)]d[m] 



I m 



I 



+ 2 



S2l 



«a(t)«*(t) 



5a(t)<5*(t) 



a2f 

(S*(t))2 



diammt)] 



(37) 



Since 



(it 



52 

ai2 



5a? ~ J * [9a2 "2 dtV^aSdJ ^ 2 0(^2 

3 d ( d'^L \ 



2 dt \ d<i>d<b ) ~^ 2dt'^ 



do? 



(38) 



and the equation of motion for a (i) implies that f^ — ^ (ff-) = Z^„>i fn{o-,d,^,^){a(t) — aci{t))'^, 



d fdL d fdL\\ ^ dUa,d,^,^) ^ 

^ ^ ' ' n>l 



5$ 



-fl^ - = E ^/-("'"'^'^) (a(t)-a,Kt))- 

n>l 



(39) 
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both vanish when a(t) = aeKt). Similarly, £ (f§ - | (§)) = ^ (f - | (§)) = 
when = ^cl{t). Prom the equation (4), 



S^L _ 3 d_ / a^L \ 
'da? 2di \dada J 



1 (i2 /a^L 



2dt^ \ da 



1 ^2 fd^L 



3a62-6ay($) + 12a + ^ ^-3e-*Md 

/73(F) 



di3 



+ e-*$3d 



-a^y"($) + 2— $d(d^ + 3ir) 
+ 3^$d(d2 + m) - 9^d(d2 + K) 



9i 
9 

+ -(200^ + 3a2d). 



51 



(40) 

and substituting the approximate solutions to the equations of motion for a{t) and <E>(i) [10], based 
on the heterotic string potential, * 



a{t) 

m 



aoe 



In 



At 

27 h'^e~^'^" 

4^ ?.3^4 cos/i(VC7(t-to)) 







fee" 



(41) 



gives 

2 \ 9a5d 



1 d2 /d^L 



d'- 



2 df^ \ dd? J t ^ CO 



3ao(l + 4A^) + 



/,2 — 3cro 

2\ I 



16 



a$2 



^_3d/^\ ld^/a^\ ^ ^a3;,2(2 + 3AV^* + ofe(3A-Vc)*)>0. 

(42) 

To second order, the probability distributions about the classical paths in minisuperspace would be 
Gaussian, and average values {a{t)), ($(t)) equal aci{t), ^ci{t)- A shift in the expectation values 
(a(i)), ($(i)) only arises at third order in the expansion of the functional integrals in equation 
(36). 
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